MINIMIZING THE GROUND STATE ENERGY OF AN ELECTRON 
IN A RANDOMLY DEFORMED LATTICE 



JEFF BAKERi, MICHAEL LOSS^ AND GUNTER STOLZ^ 



Abstract 

We provide a characterization of the spectral minimum for a random Schrddinger operator 
of the form H = —A + J^iez^ ^(-^ ~ ^ ~ ^i) L^(R°'), where the single site potential q is 
reSection symmetric, compactly supported in the unit cube centered at 0, and the displacement 
parameters Ui are restricted so that adjacent single site potentials do not overlap. In particular, 
we show that a minimizing configuration of the displacements is given by a periodic pattern 
of densest possible 2'^-clusters of single site potentials. 

The main tool to prove this is a quite general phenomenon in the spectral theory of Neumann 
problems, which we dub '^bubbles tend to the boundary." How should a given compactly 
supported potential be placed into a bounded domain so as to minimize or maximize the first 
Neumann eigenvalue of the Schrddinger operator on this domain? For square or rectangular 
domains and reflection symmetric potentials, we show that the first Neumann eigenvalue is 
minimized when the potential sits in one of the corners of the domain and is maximized when 
it sits in the center of the domain. With different methods we also show a corresponding 
result for smooth strictly convex domains. 

1. Introduction and Main Results 

1.1. The Displacement Model. The one electron model of solid state physics describes 
the behavior of a single electron moving under the presence of an exterior force generated 
by the effective potentials of a fixed configuration of nuclei in a solid. Also disregarding 
electron-electron interactions, this results in the one-electron Schrddinger operator 

H = -A + V 

in L'^{M.'^), where —A and V are the kinetic and potential energy of the particle, respectively. 

Typically one chooses a potential V that effectively models the characteristics of a particular 
solid. For example, one might use a periodic potential to model a crystal or other well ordered 
substance. While for a material containing a sufficient number of impurities, or disorder, one 
might use a random potential. In this paper we consider a potential generated by identical 
atoms or ions located at the points i + uJi, i & Z*^, i.e. 



l-UJi). 

We refer to q as the single site potential and consider real valued single site potentials q G 
which are reflection symmetric, i.e. symmetric in each variable with the remaining variables 
fixed, and compactly supported in the unit cube Aq := (— |, \Y ^"^i i-^- supp q C [— r, rY C 
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Aq, r < 1/2. We denote the collection of displacements hj u = {cUjjjg^d, where each Ui G 
[~dmax,dmaxY- Finally we choose r + d^ax = ^, which insures that adjacent single site 
potentials in the sum above do not overlap. 
For any possible collection of displacements u, 

H{uj):=-A + V^ (1) 

with domain H'^{M.'^), the second order Sobolev space, defines a self adjoint operator in L^(R'^). 
We will refer to the family H{uj) as the displacement model. As V is uniformly bounded with 
respect to lo, the spectrum of H{u), a{H{uj)), is uniformly bounded from below. 
The question we will address is the following: How can one characterize 

Eo:=mimfa{H{uj)), (2) 

i.e. the infimum of the ground state energy of H{u!) for all possible nuclear configurations tu? 
In particular, is there a minimizing configuration co""*"^ such that 

inf(T(/7(cj™")) = Eo, (3) 

and how does it look like? 

Our main result is the answer to this question: 

Theorem 1.1. A minimizing configuration cj'™" for the ground state energy of H{uj) in the 
sense of ([3]) is given by 

cur = ii-lY'dmax, i-lY'dmax) (4) 

for all i = {ii, . . . , id) G Z'^. 

The energy minimizing potential V^min is 2-periodic in each coordinate and given by the 
densest possible cluster of the nuclei in the period cell ( — |, |)'^, namely all single site potentials 
within the cluster move as close to the center (i, . . . , i) of the period cell as possible, see 
Figure [H Thus, by Floquet-Bloch theory [20], Eq is the lowest eigenvalue of —A + V^min 




Figure 1. The support of V^min for d = 2 and radially symmetric /. 

restricted to L^((— i, l)'^) with periodic boundary conditions. Due to the symmetry of the 
2'^-cluster this is the same as the ground state energy of —A + q{x — {dmax, ■ ■ ■ ,dmax)) on 
L^(Ao) with Neumann boundary conditions. For d = 1 Theorem 11.11 has been conjectured 
and partially proven in |17j . 
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Our interest in this result is mostly motivated by the case where the displacements uJi,i E Z"^, 
are independent, identically distributed R'^-valued random variables with [—dmax,dmax]'^ as 
the support of their common distribution /x. In this case we refer to H{uj) as the random 
displacement model, which is ergodic with respect to shifts in Z'^. Thus its spectrum is almost 
surely deterministic, i.e. there exists S C M such that 

a{H{ijj)) = S for a.e. ur, 

e.g. H. 

Compared to other prominent models of random Schrodinger operators, e.g. the Anderson 
model or Poisson model, few rigorous results are known for the random displacement model. 
This is mostly due to the fact that H[uj) does not depend monotonically (in form sense) on 
the random parameters uji. Even the structure of the almost sure spectrum S is unclear. It 
can be said that 

S= U criHiu)), (5) 

UJ^Cper 

where Cper is the set of all configurations u : Z'^ ^ supp /i which are periodic with respect 
to some sublattice of Z*^. This follows by adapting the proof of the corresponding result for 
Anderson models, given e.g. in [14J. 
A consequence of Theorem 11.11 is 

Corollary 1.2. The infimum of the almost sure spectrum of the random displacement model 
is given by 

inf E = Eo = inf 

Proof. This follows from ([5]) since u;™" g Cper and, by Theorem II. H inf a{H{u"^^"')) < 
inf a{H{uj)) for all other u G Cper- C] 

Note that, at least for sign-definite q, the answer to the same question for the Anderson or 
Poisson model is quite straightforward and found by considerations involving not much more 
than minimizing the potential energy: For the Anderson model, inf S is found by choosing all 
random couplings minimal. For the Poisson model one has inf E = if g > and inf E = — oo 
if g < 0. In fact, the latter, with few exceptions, only requires that the negative part of q 
doesn't vanish [1]. For the Anderson model with sign-indefinite q the description of E or just 
inf E causes difficulties similar to those for the random displacement model. Najar [18] has a 
result for this case in the small coupling regime, proven by perturbative arguments. 

What makes our result, as well as the techniques in its proof, rather interesting is that min- 
imizing the spectrum in the displacement model requires an understanding of the interaction 
between kinetic and potential energy. Physically, one can understand our result best for the 
case of negative potential wells q. In this case the formation of clusters of 2^^ sites allows for 
states with low potential energy without sacrificing much kinetic energy. But we stress that 
Theorem 1 1 . 1 1 and Corollary 11.21 hold without any sign-restriction on q. 

For the multi-dimensional random displacement model it is not yet known if the spectrum 
is localized, in the sense of being pure point, near the bottom of the spectrum. This is in 
contrast to the situation for Anderson and Poisson models. For the Anderson model this is a 
long standing result, with the hardest case of Bernoulli distributed random couplings recently 
settled in [2]. The new type of multi scale analysis introduced in [2] was now also used to 
prove the corresponding fact for the Poisson model in arbitrary dimension [HI [ID]- For the 
one-dimensional displacement model, localization at all energies was proven in ^ and, with 
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different methods and under more general assumptions, in [5]. The only available result on 
localization for the multi-dimensional displacement model is Klopp's work [13], establishing 
the existence of a localized region for the semiclassical version —h'^A + of ([T]) if /i is 
sufficiently small. 

Theorem 11.11 should serve as a first step towards understanding the spectral type of —A + V^j 
near inf S by identifying the periodic configuration in which inf S is attained. An important 
further step towards localization would be to quantify probabilistically how many other config- 
urations have ground states close to inf E, that is, to prove smallness of the integrated density 
of states (IDS) near inf E (or a related finite volume property). To this end it is interesting 
to note that for d = 1 the configuration given by (jl]) (in this case "dimerization" ) is only one 
of many minimizing periodic configurations. This will have interesting consequences for the 
IDS. In particular, one may not find the Lifshitz tail behavior familiar from Anderson and 
Poisson models and the exact asymptotics may depend strongly on the displacement distri- 
bution fi. However, we believe that, under suitable assumptions on q suggested by our proof 
of Theorem II. H in d > 2 the configuration (jl]) is the unique periodic minimizer. We plan to 
investigate this further in a separate work. 

1.2. Bubbles tend to the boundary. Theorem 11.11 amounts to optimizing the infinitely 
many parameters tu,, i G Z'^, with respect to minimizing the spectrum. Surprisingly, as will 
be shown in Section [21 its proof can be reduced to the following spectral optimization result 
in just one parameter. 

Theorem 1.3. Let q be as above, i.e. hounded, reflection symmetric, and supported in [— r, r]'^ 
for some r < \. Let d^ax = \ -r and for a e [-dmax, dmaxY Hj^^{a) := -A + q{x - a) 
in L^(Ao) with Neumann boundary conditions on dAo and denote the ground state energy of 
iJ^(a) by Eo{a). Then we have the following alternative: Either 

(i) Eo{a) is strictly maximized at a = and strictly minimized in the 2'^ corners 

maxi • • • ; i'J^max) of [ dmaxi dmax\ 

or 

(ii) EQ{a) is identically zero. In this case the corresponding eigenfunction is constant 
outside of the support of q. 

In fact, in will show that in case (i) the function -E'o(a) is partially strictly decreasing away 
from the origin, i.e. that whenever all but one of the variables (ai,...,ad) are fixed, then 
E{ai, . . . , ad) is strictly decreasing for the remaining variable in [0, 1/2 — r] and, by symmetry, 
strictly increasing in [— 1/2 + r, 0]. 

A sufficient, but far from necessary condition for case (i) to hold is that q has fixed sign 
and does not vanish identically, as in this case Eo[a) never vanishes. Case (ii) happens if the 
Neumann problem for —A + g on the support of q has lowest eigenvalue 0. Non- vanishing q 
with this property are easily constructed. 

We find Theorem 11.31 quite interesting for its own sake, independent of its application to 
prove Theorem ll.li It is a prototype of what seems to be a very general phenomenon appearing 
for Neumann problems on bounded domains, namely that "bubbles tend to the boundary". 
To this end, we have the following result for general strictly convex smooth domains and 
smooth potentials, proven in Section O with a method very different from the one we use in 
Section m to prove Theorem 11.31 
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Consider an open, bounded domain D cM.'^ with smooth boundary. We shall assume that 
D is strictly convex. Let q{x) be any bounded smooth potential whose support is a subset of 
D. For a G M*^ let qa{x) := q{x — a). In D consider the Schrodinger operator 



with Neumann boundary conditions on dD (where restriction of Qa to D is implied). We 
denote its ground state energy by -E'o(a)- As shown in Lemma 12.11 of Section [2], Eo{a) is 
continuous in a. 

Denote by G C M'' the collection of vectors a such that Qa has its support also in D. Note 
that G is an open set. 

Theorem 1.4 (Strong minimum principle for Eq). If Eo{aQ) = infagc £'o(a) for some 
ao G G, then Eq^o) is identically zero. In this case the wave function is constant outside 
the support of the potential. In other words if EQ{a) does not vanish identically in G, then 
Eo{ao) > inf^eG ^o(a) for all UqEG. 

The continuous function Eq : G M. must assume its minimum. By Theorem II. 4^ if Eq^q) 
does not vanish identically, the minimum must be assumed on dG. In the same situation 
Theorem 11.31 gives the more precise result that the minimum is assumed in the corners of 



For radially symmetric q and various types of domains D, the question of minimizing the 
first Dirichlet-eigemralVie of —A + on D is well studied, see [12] for disks and regular 
polygons, or pTj for a more general class of domains which have a certain reflection property 
with respect to the symmetry axes of the potential. 

Common to all results for Dirichlet problems is that the maximizing and minimizing posi- 
tions depend on the sign of the potential. For an obstacle g > (or "g = cx)" , meaning a hole 
in the domain marked by an additional Dirichlet boundary condition) the maximizing position 
is in the "center" of the domain, while the first eigenvalue is minimized when the obstacle is 
in contact with the boundary. The reverse is true for the case of a well g < 0. As pointed out 
in [11], this is most easily understood, if not proven, by a perturbative argument: Consider 
—A + Ag on Z) with Dirchlet boundary condition. Its lowest eigenvalue Eq{X) satisfies 



where Uq is the ground state eigenfunction of the Dirichlet Laplacian. Thus Eq{\) changes the 
most (least), if q is placed where Uq is largest (smallest), which is near the center (boundary) 
of D. The sign of q determines the sign of E'^ and thus reverses the role of maximizer and 
minimizer. 

This motivation through first order perturbation theory fails for the Neumann problem. In 
this case ([6]) still applies, but the ground state of the Neumann Laplacian is constant and thus 
Eq[0) is independent of the placement of q. This explains why the Neumann version of the 
problem is more subtle than the Dirichlet problem (roughly by one order of perturbation 
theory). Consequently, the methods used in [12] and [11] do not extend to give similar 
results for the Neumann case. An exception is a remark in [12] concerning infinite spherical 
obstacles in spherical domains. The only other work on the Neumann case, which we found in 
the literature, is [TB], which gives perturbative and numerical results concerning the optimal 
configurations of small Dirichlet holes in planar domains for maximizing the first Neumann 
eigenvalue. 




A + g, 



'a 





(6) 
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We indeed use a second order perturbation theory formula as the starting point of the proof 
of Theorem 11.41 see fl2^ below. Our proof of Theorem 11.31 doesn't use perturbation theory 
(but Floquet-Bloch theory, and the variational characterization of ground states, convexity of 
the kinetic energy, and unique continuation of harmonic functions). Still, the result may be 
motivated by second order perturbation theory: 

In d = 2 (for simplicity) consider the Neumann problem —A + Xq on L^((— i, ^y,dxdy). 
The lowest eigenvalue Eo{X) satisfies the second order perturbation formula 

s;'(o) = -25:^1^. m 

where Ek and Uk are the higher eigenvalues and eigenfunctions of the Neumann-Laplacian, 
see Section [2l3l Considering only the leading term of ([7]), corresponding to Ei = E2 = n'^ , we 
get 

q{x, y) sin^TTx) dx dy^ ^ ^) '^^ '^^ 

which is negative, independent of the sign of q. If g = g^, with go reflection symmetric and 
of fixed sign, then both integrals are zero for a = 0, and both integrals become maximal (in 
absolute value) if a is located near one of the four corners (±|, ±|). Again, this is independent 
of the sign of g. 

2. Preliminaries 

Throughout this section i7 C M'^ will be open and bounded. The Neumann Laplacian — 
on Q is the unique self adjoint operator whose quadratic form is 

\Vfix)\'dx 

for / in the domain H^{Q), the first order Sobolev space. 

2.1. Continuity of Eigenvalues. Assume that Q satisfies the i7^-extension property, i.e. 
there exists a bounded operator E : H^{Q) —>■ H'^{M.'^) such that {Ef){x) = f{x) for all 
/ e H^{Q) and almost every x & Q. Note that a sufficient condition for this is that fl has 
Lipschitz boundary, e.g. Theorem V.4.12 of [6J. 

Let g G L°°(]R'^) be real-valued and define qaix) = q{x — a) for all a G M.^. Let (a) = 
-A^ + qa. 

Lemma 2.1. (a) has purely discrete spectrum consisting of eigenvalues EQ{a) < Ei{a) < 
• • • < En{a), counted with multiplicity, where all functions En are continuous in a. 

Proof. Fix C > ||g||oo- From the extension property of Q and boundedness of g it follows that 
(-A^ + C)-i and {Hi{{a) + C)-^ are compact, e.g. Theorem V.4.13 of [6J. Thus ifj^(a) has 
purely discrete spectrum -E'o(o) < Ei{a) < . . .. It remains to show that {H^ {a)+C)~^ is norm- 
continuous in a. Continuity of the eigenvalues of {H^ (a) + C)~^, and thus the eigenvalues of 
(a), then follows from the min-max-characterization of eigenvalues. 
Without restriction, consider continuity at a = 0. For any p G {d, 00) with p > 2 one has, 
e.g. Theorem 4.1 of [21], 

||Xn(ga-g)(-A + C)-V2|| < ^2n)-''/^Xniqa-q)\m-\' + C)-'/\ 

asa^O. (8) 
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As E ■ H\n) H\R'^) is bounded, it follows that (-A + C)^/^E{-A^ + C)-^/^ jg 
bounded from L^(r2) to L'^(R.'^). Combined with ([8]) this yields 

\\Xniqa-q)i-A^ + C)-'/'\\^0 as a 0. (9) 

Norm-continuity of {Hq (a) + C)^^ at a = now follows from IQ, boundedness of (— + 
Cy/'^iH^iO) + and the resolvent identity 

«(a) + C)-' - « (0) + C)-' = (H^ia) + Cr\n{qa - g)«(0) + C)-\ 

□ 

2.2. Positivity and non-degeneracy of the ground state. We will frequently use that 
for the domains considered by us and bounded potentials q the ground state energy Eq of 
H = — + g is non-degenerate and that the corresponding eigenfunction can be chosen 
strictly positive. This generally holds if, in addition to the assumptions from Section 12.11 Q 
is connected. The latter guarantees that the ground state energy of — A^ is non-degenerate 
(— Aq = implies that | Vv?P dx = 0, i.e. Vip = and thus ip constant by connectedness). 
Non-degeneracy and positivity of the ground state of H follows from the general theory of 
positivity preserving operators provided in Section XIII. 12 and the following Appendix 1 of 



2.3. Perturbation formulas. For completeness, let us briefly recall the derivation of the 
eigenvalue perturbation formulas which we use in our arguments. Most significantly, this will 
be the first and second order perturbation formulas fl3T|) and fl25l) with respect to displacements 
of the potential in Section O However, they follow in the same way as the corresponding 
formulas for coupling constant dependence, e.g. ([Tj), so we will focus on the latter. 

Let fl and q satisfy the assumptions of the previous two subsections, H{X) = — A^ + Xq, 
Ek = Ek{X) its eigenvalues ordered hj Eq < Ei < E2 < . . . and Uk = Uk{-, A) corresponding 
real normalized eigenfunctions. Then 

= (10) 

and 

^.(^)^_2yK^_ (11) 

The formulas (13T1) and (!25|) below follow with the same argument, using smoothness of q and 
differentiating separately with respect to each component of a (the extra term (mq, (Ago)Mo) 
in (I25II does not appear in (|TT1) as (?^(Ag) = 0). 

The equation ( JTOl) is the classical Feynman-Hellmann formula, derived by using non-de- 
generacy of Eq (and thus analyticity of Eq and uo in A) and the fact that (uq, dxUo) = 0. 
Differentiating ffTOl) and using completeness of the Uk we get 

E'^{X) = 2{dxUo,quo) 

= 2^{quo,Uk){dxUo,Uk), (12) 

fc>0 

noticing that the A; = term vanishes. Differentiating the eigenvalue equation —Auq + Xquo = 
EqUq yields, for every A; > 0, 

quo = {Eo - Ek)dxUo + (mq, quo)uo - (-A + Xq - Ek)dxUo, 
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and thus 



dxuo 



quo - {uq, quo)uo + (-A + Ag - Ek)dxUo 



Eq — Ek 

After noting that ((—A + Ag — Ek)dxUQ,Uk) = 0, (ITT]) follows from inserting (fT3|) into (fT2|) . 

3. Theorem 11.31 implies Theorem 11.11 



(13) 



Theorem 11.31 says that infai?o(a) = -Eo(fl™"), where a™" corresponds to one of the 2*^ 
corners of the the cube [— rfmax, dmaxf, say a™" := {dmax, • • • , dmax)- Once we know this, then 
the central ideas of the proof of Theorem 11.11 are (i) Neumann bracketing to go from H{uj) to 
operators of the type H^^ (a) and (ii) extending the ground state of the minimizer H^^ (a*"*") 
to by repeated reflection. 

Proof, (of Theorem I l.ip For any given configuration a;, the restriction of Hiuj) to the unit cube 
centered at z G Z*^ with Neumann boundary conditions is unitarily equivalent (via translation 
by i) to Hj!^^{uJi), defined as in Theorem 11.31 Thus, by Neumann bracketing and Theorem II. 3 [ 

\iezo J 
> inf {Eo{a) : a G [-dmax, dmaxY} 
= Eo(a™"). 

This holds for arbitrary configurations uj and thus, hj , Eq > EqIq"^^^). 
Now consider cj™*"- = (ct;™")^^^^ as given by (jlj). The corresponding potential 

is 2-periodic in Xi for each i. By Floquet-Bloch theory [20] the bottom of the spectrum of 
if (a;™") = —A + V^min is given by the smallest eigenvalue Eq'^^ of its restriction to Aq := 
(— |, with periodic boundary conditions, see Figure [21 













Figure 2. The period cell of V^^ir, in d = 2. 

On Aq the potential V'^"^'" is symmetric with respect to all hyperplanes Xi = 1/2, i = 
1, . . . ,d. Thus Eq'^^ coincides with the smallest eigenvalue of the Neumann problem on Aq. 
Again by symmetry of the potential, the latter coincides with the smallest eigenvalue of the 
Neumann problem on Aq. As cj™" = {dmax, ■ ■ ■ , dmax) = a™", this eigenvalue is £'o(a™"). In 
summary we have shown that 

Eo < inf a(/i(cu™")) = Eo(a™") < ^o- 
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Thus Eq = inf (T(if(a;™")), which proves Theorem II. 1[ □ 

4. Proof of Theorem 11.31 

This entire section is devoted to prove Theorem 11.31 Thus we work under the assumptions 
that q G is real-valued, non-vanishing, reflection symmetric and supported in [— r, r]*^, 
r < 1/2. 

Suppose that alternative (ii) of Theorem 11.31 is false. We will show that this implies that 
alternative (i) must hold. 

We begin by fixing all of the components of the displacement parameter except for one, 
which may be chosen to be the first, and consider the lowest Neumann eigenvalue as a function 
only of the first coordinate, i.e. Eq^g) := i?o(a, ^2, • • • , ^d)- We note that Eq^q) depends 
continuously on a (see Lemma 12. ip and that by symmetry we have 

Eo{-a) = Eo{a). 

For this reason we will restrict ourselves to the case a G [0, | — r] and show that 

Eo{ai) > Eo{a2) for < ai < a2 < 1 - r. (14) 

As the same holds for Eq as a function of each other coordinate, we conclude from this that 
Eq has a strict maximum at the origin and strict minima at the corners {±dmax, ■ ■ ■ , ^dmax), 
i.e. we are in the situation of alternative (i). 

As (a2,...,arf) will be kept fixed, we will use the (slightly sloppy) abbreviation q{x) : = 
q{xi, X2 — a2, . . . ,Xd — ad) for the rest of the section. For a scalar a and ei = (1, 0, . . . , 0) we 
then write qa{x) = q{x — aei), a notation to be used also for functions other than q. By dn 
we will denote the exterior normal derivative on the boundary of a given domain. 

Lemma 4.1. Eq{0) > Eq^o) for every a G [0, | — r]. 

Proof. Define the tube L := {x : \xi\ < = 2,. . . ,d} and construct a periodic extension, 
qper{x), of the potential q{x) on L by 

qper{x) = y^^qi{x). 

We consider the Neumann problem on L for —A -|- qper{x), see Figure O Let ip denote 



N 



® 


® 




® 



N 



Figure 3. Periodic extension, qper{x), of g(a;) 

the normalized Neumann ground state of —A -|- q{x) on the unit cube Aq. Construct a new 
function \1/ on the tube L by periodically extending ip on all of L. Symmetry of the potential 
implies that is a smooth solution of 



- A^(X) + qper{x)ij{x) = Eq{0)^{x) 



(15) 
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on all of L. Now multiply (|T5|) by '^{x). Then over any unit cell C in L, periodicity of ^ 
implies we may integrate by parts without creating boundary terms. In particular, this holds 
for the unit cell := Aq — aei, yielding 

^o(O) = / {\V^{x)\^ + qper{x)^>\x))dx. (16) 

JKa 

Shifting to the right by a does not affect the result of equation (fT6|) . i.e. 

^o(O) = / {\V^a{x)\'' + qper{x-aei)^l{x))dx 
Jao 

= [ {\V^,{x)\' + qa{x)<ill{x))dx. (17) 

While "^aix) does not satisfy Neumann boundary conditions on Aq, it is still in the form 
domain H^{Ao) of the Neumann operator —A + qa{x) on Aq. Therefore, minimizing the right 
hand side of equation f[T7l) over all normalized functions in H^{Aq) it is clear Eq{0) > Eq^o). 

To show that indeed -Eo(O) is strictly greater than EQ{a) it suffices to show that, when 
restricted to Aq, is not equal to a multiple of the Neumann ground state eigenfunction 
corresponding to the potential qa on Aq. Suppose, for contradiction, that was such a 
multiple. Then by construction the box B = (— i, — ^ + a) x ( — i, is disjoint from the 
support of the potential and satisfies the equation 

-A-^a = Eo{0)-^a (18) 

with Neumann conditions on the boundary of B. As \l'a > 0, it is the ground state of the 
Neumann problem on B. Thus Eo{0) = and niust be constant on B. This entails that 
is harmonic outside the support of the potential, and since it is constant on an open subset, by 
unique continuation of harmonic functions it must be constant everywhere outside the support 
of the potential. However this implies that alternative (ii) must hold, a contradiction. Thus 
Eo{0) > Eoia). □ 

Lemma 4.2. For any positive integer n and a G (0, | — r], Eo{na) > Eq^^u + l)a) so long as 
{n + l)a is less than or equal to \ — r. 

Proof. To keep notations simple, we first show this for n = 1. Again consider the tube 
L := {x : \xi\ < = 2,...,d}. Fix a with < 2a < | — r and consider a 2-periodic 
extension, Wper, of the potential 

w{x) := qa{x) + q-a+i{x) (19) 

on L given by 

Wper{x) := y^W2i{x). 

As before we consider the Neumann problem on L for —A + w.p(.r[x\ see Figure HI and let ^/^ 
denote the Neumann ground state of — A + ga(x) on the unit cube Aq, normalized to = |. 
Construct a new function \E' on all of L by 2-periodically extending ip{xi,X2, ■ ■ ■ , Xd) +ip{—xi + 
l,a;2, . . . ,Xrf). 

Symmetry and periodicity of the potential implies that \E' is a smooth solution of 

-A^(x) + Wper{x)'^{x) = Eo(a)^(x) 
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Figure 4. Periodic extension, Wper{x), of w{x) := Wa{x) + w-a+i{x) 



on all of L. Now we proceed in analogy to f|T6l) and f|T71) . this time considering cells of length 
2 instead of unit cells and again shifting by a to the right. We find 



Eo{a) 



-a<xi<i 



-a} 



|V^(a;)p + Wper{x)-^{xf) dx 



{x'eL;-|<xi<|} 



(20) 



As above we conclude that Eo{a) is not smaller than the first Neumann eigenvalue E of 
—A + g2a + ?i on {a; G L : — I < Xi < |}. A further application of the argument in 
the last paragraph of the proof of Lemma 14.11 shows that \E'a restricted to the unit 2 cell 
{x e 1/ : — i < Xi < |} is not a multiple of the ground state eigenfunction corresponding to 
E. Thus Eola) is strictly greater than E. 

Imposing an additional Neumann condition at Xi = | can not increase the lowest eigenvalue, 
thus 

Eo{a) > E> min{Eo(2a), ^o(O)} = Eo{2a) 

by Lemma 14. H which concludes the proof for n = 1. The crucial idea which allowed us to 
reduce the n = 1 claim to Lemma 14.11 was that the term q-a+i in 0191) was shifted back into 
the center of the cube {x G L : | < xi < |} in (|20l) . The same mechanism can now be used 
to inductively prove the claim for all n. □ 

We can now readily complete the proof of f|T^ : As a above was arbitrary. Lemma 14.21 
implies that Eo^a) is strictly decreasing on the set of all dyadic numbers in [0, | — t]. As this 
set is dense in [0, | — t] and E^^a) is continuous, then Eq^o) is strictly decreasing on all of 
[0, 1 — t]. Therefore assuming (ii) is false, (i) must be true. 

5. Proof of Theorem 11.41 



In Theorem 11.41 it is assumed that the domain and the potential are smooth. Thus we have 
by elliptic regularity that C°°{D) is a form core for the Neumann operator H^{a) = —A^ + qa- 
Moreover, the eigenfunctions are all in C°°{D) and have normal derivative zero on dD. 

We call the eigenvalues Ek{a), ordered and accounting for multiplicity, and the normalized 
eigenfunctions Ufc(x; a), k = 0,1, . . .. We choose the ground state Uq strictly positive and all 
other Uk real. They form an orthonormal basis in L'^[D). Thus, for any function / G H^[D) 
we have that 



k=0 



D 



0V/(x)p + g,(x)/(x)2] dx 



(21) 



where (■,■) denotes the inner product on L'^[D). Note that / is not in the domain of the 
operator, just in the form domain. See [20] and [7]. 
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Let N denote the outward normal vector field on dD. It will be convenient to use that 
can be extended to a smooth vector field in a neighborhood of dD. To see this, first work in a 
neighborhood of a fixed point of the surface. Without loss (i.e. up to a rigid motion) we can 
choose this point to be the origin and the surface to be given by 



We get a global extension of iV to a neighborhood of dD by using compactness of dD and a 
standard partition of unity argument. 

In this neighborhood we define the matrix- valued smooth vector field K = (Kij) = {diNj). 
The restriction of K to dD is the curvature matrix of the surface. Indeed, in the local 
coordinates used above, we have 



We have assumed that D is strictly convex. This means that the Hessian of / is negative 
definite and thus at every p e dD the restriction of K{p) to the tangent plane at p is positive 
definite. 

The following identity for Eq^o) is the main technical ingredient into our proof. Here AEq 
and V-Eq refer to the a-derivatives of Eq. Otherwise, all symbols such as 9j, V, A denote 
derivatives with respect to the spatial variable. B{-,-) is the bilinear form 



Lemma 5.1 (Second order perturbation theory). The ground state energy satisfies the 
equation 



Xd = f{xi, . . .,Xd-l) 

in a vicinity of the origin, where / is smooth, 

/(0,...,0) = and V/(0,...,0) = 0. 
Thus at a point p = (xi, . . . , Xd-i, f{xi, . . . , Xd-i)) near one has 



(22) 



(23) 





B{u,v) = {u,Av) — {Au,v). 




2 




(24) 



Proof. We start with the second order perturbation theory formula 




(25) 



see Section 12.31 Differentiating the eigenvalue equation yields 

(Vga)Mo = ^oVmo - (-A + qa)Vuo 



and therefore 



(wfc, (Vga)Mo) 



{Ek - Eo) {uk, Vmo) + B{uk, Vmo) • 



(26) 
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Hence 

AEo = {uo,{Aqa)uo) ~2j2[iEk- Eo){Vuo,Ukf -2{Vuo,Uk) ■ B{uk,Vuo)] 



Ek — En 



which, using (!26|) once more, can be rewritten as 

AEo = {uq, (Ag„)no) + 2 Vuq) + (n^, (Vg„)?io)] • Vuq) 



-2E 



-fi't — E(] 



(uo, (Aga)tto) + 2 Vuo) + {uk, (Vga)Mo)] ■ (wfc, Vuo) 

A: 

- 2[5(mo, Vmo) + {uq, (yqa)uo)] ■ {Uq, VUq) 
Y,.B{uk,diUQY 



-2E 



El. — En 

= (mo, (Aga)Mo) + 2(Vmo, (Vga)Mo) + 2 ^ Vmq) ■ (wfc, Vmq) 
- 2[5(mo, Vmo) + (mo, (Vga)Mo)] ■ (^io, Vmo) 

^ Efc - 

where finally the completeness relation of the Uk was used. It is clear that 

(mo, (Aga)Mo) + 2(Vmo, (Vga)Mo) = 
and using again (126|) with A; = we can simplify further and get 

AEq = 2 ^ Vuq) ■ {uk, Vmo) - 4(uo, (Vga)Mo) ■ (^io, Vmo) 

k 

^iB{uk,diUof 



-2E 



— En 

Recall that B{u, v) = {u, Av) — {Av, u) and hence 

B{uk, Vmo) ■ (wfc, Vmo) = ^[(wfc, AVmo) - {Auk, Vmo)] ■ (wfc, Vmq) • 

Since uq G C°°{D) we know that the vector AVmq has square integrable components 
hence 

^(ufc, AVmo) • (ufc, Vmo) = (Vmo, AVmo) • 

k 

The second term of (1281) we write as 



^([-A + qa]Uk, VMo) ■ (Wfc, Vmo) - ^(gaMfc, Vmo) ■ (Wfc, Vmo) . 
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The second sum equals ^j{djUo, qadjUo) while the first sum is 

j k j 

since djUo is in the form domain. Collecting terms we find from ( |28l) that 



V / {djUo)N-V{d,uo)dS, (29) 



where Green's identity was used. On an open neighborhood of dD we have 

J2id,uo)N-V{djUo) = J2{d,uo)d,N-Vuo-J2idjUo){djN,){d,uo) 

j j hi 

= Vuo ■ V{N ■ Vmo) - Vmo ■ KVuo , (30) 

where K is the curvature matrix defined above. Using that ■ Vuq = on dD one has that 
the first term is VjMo ■ Vt(A^ ■ Vuq) for points in dD, where denotes the component of 
the gradient in the tangential directions. However, ■ Vuq = and smoothness of dD also 
implies Vt{N ■ Vuq) = and thus Vmq ■ V(A^ ■ Vuq) = on dD. Thus and yield 



B{uk, Vuq) ■ {uk, Vmo) = - / Vmo ■ KVuo dS. 

^ JdD 



After substituting this and the first order perturbation formula 

V^o = -{uo, (VgjMo) (31) 

into ([27]) we arrive at □ 

Lemma 5.2. Assume that q is a smooth potential with compact support in D. Assume that 
the Neumann ground state Uo{x;a) for some fixed E G is constant on dD. Then there 
exists an open neighborhood of Uq where £'o(ao) ^ E^la). 

Proof. By shifting coordinates we assume that = 0. We shall proceed by a trial function 
argument. Consider the problem with the shifted potential qa{x) := q{x — a). Denote by 
Da = D + a the shifted domain, i.e., the function Ua{x) := Uo{x — a; 0) solves —Aua + qaUa = 
Eo{0)ua, with a Neumann boundary condition on dDa as well as being constant on dDa- We 
shall construct a trial function in the following fashion. In the intersection of Da with D 
we set (p = Ua and in D\Da we set to be a constant which equals the boundary value of Ua- 
Note that G H^{D). By the variational principle 

^^^^^ ^ |V0P + qa<P'dx 
~ Jj-^ (jP'dx 

For a E G the right side equals 

iDnPa l^^^l' + (l-^'ldx ^ \Vua\^ + qguldx ^ ^^^Q^ Ip^ldx ^^^^ 
(p'^dx ~ (p'^dx (p'^dx 

In the case Eq{0) = 0, this implies EQ{a) < 0, as was to be shown. 
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< 1 



> 1 



Next we claim that for a sufficiently small 

Id u^dx 
Jj~i (p'^dx 

if Eo{0) > and 

Id u^dx 
(p'^dx 

if Eq{0) < 0. This yields the lemma for the remaining cases. If we denote by c the boundary 
value of Mo we find that the claim follows once we show that in a vicinity of the boundary, 
Uq < c for -Eo(O) > and Uq > c for Eq{0) < 0. To see this, fix a point on the boundary, 
call it the origin and use the local coordinates (122!) and (!23l) above. The normal vector at is 
(0, . . . , 0, 1) and hence the normal derivative equals ddUo{0) = 0. Further, since Uq is constant 
on the boundary we find by differentiating uo{xi, . . . , Xd-i, f{xi, . . . , Xd-i)) = c that 

didjUoiO) = -ddUo{0)didjf{0) = 

for i, j = 1, . . . d — 1. Hence 

ddddUoiO) = Amo(O) = -Eo(0)mo(0) 

which is negative for Eq{0) > and positive for Eo{0) < 0. As this holds at all points of the 
boundary, we get the required property of Uq in a vicinity of the boundary. □ 

We remark that in the case Eo^ao) ^ 0, the above proof actually gives the strict inequality 
-E'o(ao) > E^ia) for a close to Oq, i.e. Ei^{a^ is a strict local maximum. This is the case, for 
example, if q is sign-definite, and allows for a shorter argument in the following completion of 
the proof of Theorem II. 4[ 



Proof of Theorem \1.4\ Assume that Eq attains its minimum value in G say at the point qq. 



This entails that V-Eo(cto) = and AEo{aQ) > 0. Using the Lemma [5. II we find that 

AEM = -2[ Vu,. KVu,dS -2Y, ^^l^''''^f'^' ■ 
JoD Ek-Eo 

The right side is non-positive, since D is convex. It cannot be strictly negative, since that 
would contradict the assumption that Eo^a) has a local minimum at oq. 

Thus the right side must vanish. Since D is strictly convex, the ffist term can vanish only 
if Uq is constant on the boundary (recall that K is positive definite on the tangent space 
at each point of dD and that Vuq is a tangent vector). However, Lemma 15.21 shows that 
-E'o(ao) > Eo{a) in a neighborhood of oq. Thus Eo{a) = -E'o(ao) in this neighborhood. Since 
the set where Eo{ao) = Eo{a) is closed and open in G (the above argument applies to every 
a with i?o(ao) = Eo{a)) the function EQ{a) must be constant. Assuming now that E^la) is 
constant, Uq is constant on the boundary and moreover, there must be equahty in (152]) . This 
means that 

iVual'^dx = . 



'Da\D 

Hence Ua must be constant in Da\D and therefore, for small a, = —Aua = Eo{a)ua there. 
Thus, Eo^a) = 0. Since the support of qa is a subset of D (as a G G), Ua is constant in the 
non-empty open set Da\D, which is disjoint from the support of qa- Since it is harmonic 
outside the support of qa it must be constant there too. It follows that uq is constant outside 
the support of g. □ 
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6. Discussion: Extensions and open problems 

We conclude with some remarks about possible generalizations and open problems related 
to our main results. 

(i) Theorem 11.31 and its proof immediately generalize to the Neumann problem on an ar- 
bitrary rectangular box {x : \xi\ < = 1,. . . ,d} rather than the unit cube Aq, which we 
chose to keep notations simple. This also gives a corresponding version of Theorem II. ![ where 
7j'^ is replaced by an arbitrary rectangular lattice. 

(ii) In Theorems 11.11 and 11.31 we also may replace the obstacle q by an reflection symmetric 
hole with Dirichlet boundary conditions, often interpreted as an infinite barrier. More pre- 
cisely, let C C {x : \xi\ < r,i = 1, . . . ,d} he closed and reflection symmetric. Let Ca = C + a 
and := —A on A \ Cq with Neumann conditions on dA and Dirichlet conditions on dCa- 
Then Eq^o) is minimized when the hole is in a corner of Aq and inf a{H{u)) is minimized for a 
periodic configuration of 2'^-clusters of holes. While we expect that Theorem 11.41 also extends 
to this situation, at least for holes with smooth boundary, our proof does not extend directly. 

(iii) Theorem 11.41 covers the situation of a radially symmetric potential (or Dirichlet hole) 
placed in a spherical domain, where all placements of q which touch the boundary are equiv- 
alent minimizing positions. However, in this case the methods of Lemma 15.21 may be used to 
also show that the maximal position occurs when the potential is centered in the domain. 

To see this, suppose D is a spherical domain centered at and let g be a radially symmetric 
potential compactly supported in D, also centered at 0. Radial symmetry of the domain and 
potential then imply the ground state eigenfunction, uq, corresponding to -E'o(O) is radially 
symmetric. Thus it satisfies the conditions of Lemma 15.21 If Eo{0) > one may show by 
reducing the problem to one dimension, or by using maximum principles, that outside the 
support of the potential uq < c, where c denotes the value of uo on the boundary of D. 
Similarly if Eo{0) < 0, uq > c outside the support of the potential. One may conclude using 
the arguments of Lemma 15.21 that Eo{0) > Eo{a) for every a ^ G. This then leads to the 
following alternative for the case of a spherical domain centered at with radially symmetric 
potential: either Eo{a) = 0, in which case the corresponding eigenfunction is constant outside 
the support of the potential or -E'o(a) assumes a strict maximum at -E'o(O) and strict minima 
when the support of the potential touches the boundary. 

(iv) Theorems 11.31 and 11.41 are proven with very different methods and apply to mutually 
exclusive classes of domains (rectangular boxes vs strictly convex domains). It would be 
desirable to find a method of proof which covers both results, as this method would most 
likely also cover more general polygons and polyhedra. Particularly interesting cases would 
be equilateral triangles or hexagons, as they tile the plane and would lead to a corresponding 
extension of Theorem 11.11 

(v) One can view Theorem 11.11 as a mechanism in which the nuclei of a solid self-organize 
into a simple periodic pattern, given a density condition (exactly one site per cube). It would 
be wrong, in our opinion, to see this as a model for crystallization since the regularity of the 
pattern is to a large extend determined by the density condition. Real crystallization, however, 
cannot be explained by the interaction of one electron with nuclei alone. It is a many-body 
effect and the nuclear repulsion and, more importantly, the Pauli exclusion principle play a 
role. Further, one needs sufficiently many electrons, e.g., a half filled band. Indeed, there have 
been results in this direction in [13] for the Falicov-Kimball model, a variant of the Hubbard 
model where the nuclei are treated classically and sit on a lattice where the electrons hop. 
Crystallization was shown in [13] for the half filled band. In our model, without the density 
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condition, we expect that the nuclei would stick together. While this is an open question there 
is some evidence in this direction. For bosons it was shown in [13] that the nuclei indeed stick 
together. 

It would be interesting to consider an extension of our model, a continuous analog of the 
Falicov-Kimball model, in which one considers a finite periodic array of cubes on a torus. 
Assuming the same number of spinless fermions as the number of cubes and assuming one 
nucleus in each cube, it is not unreasonable to expect that in an energy minimizing configu- 
ration the nuclei sit at the center of each cube. This is an interesting open question. Needless 
to say that the methods in this paper have no bearing on this problem. For an overview of 
the Falicov-Kimball model the reader may consult j^. 
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